Noncollapsing solution below Tc for a randomly forced particle 
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We show that a non coUapsing solution bellow r^ can be constructed for the dynamics of randomly 
forced particle interacting with a dissipating boundary. The scaling analysis predicts a divergent 
collision rate at the boundary for the non-collapsing solution. This prediction is tested numerically. 
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Recently opposing viewpoints appeared in the litera- 
ture H, H regarding the localization properties of a 
one-dimensional particle subject to an uncorrelated ran- 
dom force interacting with a dissipating boundary. The 
equation describing the particle dynamics is 



(1) 



where r]{t) is Gaussian white noise, {r]{t)r]{t')) = 2S{t — 
t'). The dissipating boundary condition is set such that 
the particle approaching with velocity — u, (u > 0), at the 
boundary a; = is reflected with velocity ru, r < 1. The 
problem was explored further with various techniques 
^, ^ , which have obtained the same critical value for the 
dissipation parameter Tc and the persistence exponent, ft 
is an unresolved problem, the absence of the collapsing 
behavior in numerical simulations 

In this paper, we propose a solution of this paradox. 
We show that for r < Vc one can construct a constant 
mass (noncollapsing) solution starting from the collaps- 
ing one. In our approach, we use the Fokker-Planck equa- 
tion (FPE) description of the process Eq. (0). The FPE 
associated with the Langevin equation (nh is 



with the dissipating boundary condition 

P{0,-u)^r^P{0,ru), u > 0, 



(2) 



(3) 



and the initial condition P{x, u,t — 0) = S(x — xq)S{u — 
Wo)- 

InRef. §, it was shown that the the general solution 
of this problem has the following integral form 

P{x,u;xo,uo;t) = Po{x,u;xo,uo;t) 
+ Jq dti duiUiP{x,u;0,rui;t — ti)Po{0,~ui;xn,uo: 

(4) 

where PQ{x,u;xo,uo;t) is the solution of the FPE Eq. 
(0) with absorbing boundary at x = M. 



Burkhardt has shown in Ref. Q that a collapsing so- 
lution with an algebraic temporal decay can be found for 
the Eq. (|). The surviving probability Q(xo,UQ,t) = 
J dxdvP{x, u; xq, UQ,t) behaves asymptotically as 



QiO,-u,t) 
Q(0,u,t) 



2 sin 



u > 0. 



u > (5) 
(6) 



For the collapsing solution (0 > 0) we have Q(0, 0, t) = 0, 
that is, the origin of the phase space (x = 0, u = 0) is 
an absorbing point for the random particle. With this 
observation the collapsing behavior can also be obtained 
numerically. Using a discretization of Eq. (|l|) together 
with the absorbing prescription at the origin, the collaps- 
ing behavior is obtained with the persistence exponent in 
perfect agreement with theoretical result of Ref. Q] , see 
Fig. In numerical simulation the particle was absorbed 
if its velocity after the collision with the boundary was 
smaller than -y/Ai, where At is the time integration step. 
As At decreases the collapsing behavior is preserved sig- 
naling the existence of the collapsing behavior in the con- 
tinuum limit. 

On the other hand, Eqs. (^), (|6|) accept the solution 
(/> = at any value of the dissipation coefficient r case in 
which the collapse does not occur since Q is constant. We 
make the second observation that Eq. (||) is not in con- 
tradiction with a solution that has constant mass, that 
is, /p dx J_^duP{x,u;xo,uo;t) = 1. Indeed if one in- 
tegrates the over x and u the right hand side of Eq. (^) 
is the conservation law for the particle in case of the ab- 
sorbing solution. The reason that allows the existence of 
more than one solution for this problem is that the dissi- 
pating condition P(0, —u) = r^P(0, ru) does not specify 
uniquely the solution on the boundary a; = 0, but is just 
a condition that the solution must obey on the boundary 
X = 0. The continuity condition for the solution asks 
^^^that P(0, It ^ 0) or oo. 

We can construct a noncollapsing solution starting 
from the collapsing one. Let us start from the FPE with 
a source term 



d_ 

di 



d \ 



/(t)<5(x-e)^(u), (7) 
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where e > and f{t) is at the moment an arbitrary 
function of time to be determined from the conserving 
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FIG. 1: First return distribution for r = 0.1, 0.05, (rc « 
0.163). The lines have the theoretical exponent found in Ref. 



mass condition. The solution of the above equation that 
satisfies the boundary condition Eq. ^ is 



Grix, u; xq, Mo; t) = Pr{x, u; xq, uq; t) 

+ f^dtiPr{x,u-e,Q-t~h)f{ti), ^ ' 

where Pr{x^u]XQ,UQ]t) is the collapsing solution of the 
problem at given r < re- After time Laplace transform, 
we have 

G{x, u] xo,uo] s) = Prix, u; xo, uo] s) + Pr{x, u; e, 0; s)f{s). 

(9) 

We can choose /(s) such that the mass of G is constant. 
We see that the choice 



/(«) 



- dxi duiPr{xi,ui]Xo, Up; s) 
dxi duiPr{xi,ui; e, 0; s) 



(10) 



gives us the needed solution. For any e > the solution 
Eq. satisfies the boundary condition since we use Pr, 
and conserves the mass by construction. In the limit e —> 
it satisfies also the initial condition 6{t)d{x ~ xo)6{u — 
Uq). We make the observation that /(s) = for r > 
Tc as Pr{x,u]Xo,uo]t) has a constant mass. It remains 
to be shown that the limit e — > exists. For the case 
r = the asymptotic expressions for Pr{x,u;xo,uo; s) 
and J_^duPr{x,u;xo,UQ;s) were obtained in Ref. 
and one can see explicitly that the above limit exists. 

We can obtain the behavior of the collision rate at 
the origin for small e using the scaling properties of the 



solution of the FPE. Equation (g) gives 

P{x, u; Xo, Uq; s) = X^P{X^x, Xu; X^xq, Xuq; A~^s). (11) 
The collision rate i?coii is given by the following relations: 

^coii(e;s)= / du uG{0,u;xo,uq;s) 



du uPr{0, U; Xq, Uq; s) 

) 

J°^du uPr(0, u;e,0;s) 



dx du Pr{x, u; e, 0; s) 

rl 

s 



dx / du Pr{x,u] xo,uq; s] 
Jo 



Using the scaling property, Eq. (|ll[), with A = 
we have for the term depending on e: 

Qi'''^) dx du Pr{x,u;e,0;s) 

J(e2/3s) 



-1/3 



,-20/3 



E2/3g(e2/3s) 



e < 1, 



(12) 



-1+0 



where we have used that J(s) « 1— cs"^ and Q{s) : 
for s < 1. 

Equation (|l^) shows that the conserving solution has 
divergent collision rate at the origin. This implies 
that the the current density uG{0, u\ xo,yo] s) is noninte- 
grable, as was first noted in Ref. [||. 

The prediction of a divergent collision rate can be 
checked numerically. If we integrate the Langevin equa- 
tion with a finite time step At, then the particle is in- 
jected at e « (At)3/2, where e is its velocity that is very 
small. Consequently, the number of bounces at the origin 
must diverge as {At)~'^ as At 0. Indeed Fig. ^ shows 
a perfect validation of this prediction. In the same figure 
we have plotted the probability that the particle stays 
between x = 0, x = 0.01. We see that this this probabil- 
ity is independent of the integration step for each value 
of the restitution parameter r. This means that there 
is no collapsing behavior. The particle is attracted by 
the wall, performs an infinite number of collisions (in the 
limit At 0) and with probability 1 is injected back 
into the domain a; > 0. 

Now we can see that the question "what the particle 
does if we put it at a; = 0, m = 0?' is indeterminate for 
r < Tc- The general solution in this case is 

P = qGr{x, U, Xq, Uq, t) + {I ~ q)Pr{x, U, Xq, Uq, t) , 

0<g< 1, 

(13) 

and one has to specify q for the answer. 

In conclusion, we have shown that the collapsing be- 
havior can be found numerically if one notice that the 
collapsing solution has an absorbing point into the ori- 
gin, hence it must be enforced in the simulation. 
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We have constructed a noncoUapsing solution for the 
case r < re- This is possible since the absorbing bound- 
ary condition P{0, —u) = r^P{0, ru) allows for two func- 
tions as boundary condition. One goes to zero as u — > 
and the other one diverges to oo as u ^ 0. 



In terms of Brownian paths, we propose the following 

picture: for r > Tc the probability to touch the origin 
of the phase space starting from any other point is zero, 
similar to simple diffusion in two or more dimensions, 
thus the collapsing solution is forbidden. When r < Tc 
the diffusing particle touches the origin with probability 
1 and if the path is set to terminate there the collapse 
occur. If the path is not set to terminate the particle is 
sent back into the domain a; > after an infinite number 
of collision with the boundary. The weight of the paths 
leaving the origin without touching the boundary is zero 
in the continuum limit but they give a finite contribution 
because they are sampled an infinite number of times. 



FIG. 2: a) Collision rate at a; = function of integra- 
tion time step At at various values of the restitution co- 
efficient r. Wc sec that for r > rc ~ 0.163 the collision 
rate is independent of At whereas diverges like (At)^** for 
r < Vc- The lines plot the theoretical prediction with 4> 
= 0.25(r = 0), « 0.087(r = 0.1), = 0(r > rc). For r = 0.1, we 
considered also the subleading correction, b) Probability to 
find the particle in the interval (0, 0.01) function of At. The 
quantity is constant for both r < rc and r > rc- The graphs 
were displaced vertically for clarity. 
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